Many-body Landau-Zener effect at fast sweep 
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The asymptotic staying probability P in the Landau-Zener effect with interaction is analytically 
investigated at fast sweep, e = 7rA 2 /(2fin) <C 1. We have rigorously calculated the value of 7o in the 
expansion P — 1— e + e 2 /2 + e 2 Io for arbitrary couplings and relative resonance shifts of individual 
tunneling particles. The results essentially differ from those of the mean-field approximation, ft 
is shown that strong long-range interactions such as dipole-dipole interaction (DDI) generate huge 
values of 7o because flip of one particle strongly influences many others. However, in the presence 
of strong static disorder making resonance for individual particles shifted with respect to each other 
the influence of interactions is strongly reduced. In molecular magnets the main source of static 
disorder is the coupling to nuclear spins. Our calculations using the actual shape of the Feg crystal 
studied in the the Landau-Zener experiments [Wernsdorfer et al, Europhys. Lett. 50, 552 (2000)] 
yield Io that is in a good agreement with the value extracted from the experimental data. 

PACS numbers: 03.65.-w, 75.10.Jm 
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I. INTRODUCTION 



Landau-Zener (LZ) effect^ (see also Refs. CTlsl) is 
a well known quantum phenomenon of transitions at 
avoided level crossing, see Fig. ^ LZ effect was encoun- 
tered mainly in physics of atomic and molecular collisions 
(see, e.g., Refs. |6||3 and references therein). In the time- 
dependent formulation, the LZ effect can be modeled by 
a two-level system (TLS) 

H = ~W(t)a z + ^Aa x , (1) 

where a a , a = x,y, z are the Pauli matrices and 

W(t) = E,(t)-£f(t) (2) 

is the time-dependent bias of the two bare (A = 0) energy 
levels. The general state of this model and the probability 
P(t) to stay in the \\) state can be written as 

*(t)=o 1 (*)U> + o T (t)|t), P(t) = \ ai (t)\ 2 . (3) 

The initial condition is W(— oo) = — oo and P(— oo) = 1. 
If W(t) changes fast, the system does not have enough 
time for transition to the state |f) and it practically re- 
mains in the state \\) , thus P(t) = 1. In the opposite case 
of slow W(t) the system mainly remains at the lower of 
the adiabatic energy levels 



1 



B ± (t) = ±-Vr(t) + A 2 , 



(4) 



thus the asymptotic staying probability P(oo) = 0. 

The time-dependent Schrodinger equation for the 
Hamiltonian of Eq. (JJJ can be solved exactly for the linear 
sweep W(t) = vt, the result for the asymptotic staying 
probability being 



P = P(oo) = e 



7rA 2 

2hv 



(5) 
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As the Schrodinger equation (SE) for a spin 1/2 is 
mathematically equivalent to the classical dissipationless 



FIG. 1: A pair of tunnel-splitted levels vs. energy bias W(t). 
Here E\, and Ef are the bare energy levels (A = 0), whereas 
E± are the exact adiabatic energy levels of Eq. (0J . P denotes 
the probability to remain in the (bare) state \[) after crossing 
the resonance. 



Landau-Lifshitz equation (LLE), the effect can be viewed 
upon as a rotation of a classical magnetization vector^S 
Recently the LZ effect was observed on c rystals of 
molecular magnets Fe 8 in Refs. fTnUTTl (see Refs. filial for 
a recent review). This posed a new problem of the many- 
body LZ effect that can be described by the Hamiltonian 
of the transverse Ising model 

i ij 

, (6) 

where Vj are the local shifts of the resonances for individ- 
ual TLSs that can be induced, for instance, by nuclear 
spinsi* and is their coupling. The SE for such a sys- 
tem of JV > 1 coupled TLSs contains 2 N time-dependent 
coefficients of the wave function. In the language of the 
time-dependent bare energy levels, there are 2 N different 
lines that cross each other at different values of W, see 
Fig. El 
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FIG. 2: Bare energy levels vs the sweep field W for a ring 
of fO two-level systems with Vi = 0, coupled by an antifer- 
romagnetic nearest-neighbor interaction. All the lines except 
the two ones with the highest slopes that correspond to the 
states with all spins up and all spins down are strongly degen- 
erate. The energy levels for the ferromagnetic coupling can 
be obtained by the upside-down transformation of this graph. 

Eq. © can be extended by taking into account 
coupling to environmental degrees of freedom such as 
phonons, as was done for one-particle LZ effect in Refs. 

One can also consider dynamics of nuclear spins- 
that are treated in a simplified way as static disorder in 
Eq. 10. As quantum dynamics of many-body systems 
far from equilibrium and coupled to environment is very 
involved, simplified theories using postulated rate equa- 
tions (neglecting quantum-mechanical coherence) have 
been proposed for both the LZ effect and the relaxation 
out of a prepared statei 18 ' 19 i 20 i 21 i 22 All these theories 
consider the molecular field at the cite i to determine 
whether the particle i is in the vicinity of the resonance 
and thus can flip. This means that these theories are 
based on the mean-field approximation (MFA). 

One can perform the MFA for the dissipationless sys- 
tem described by Eq. © by considering a single particle 
undergoing a LZ transition in the effective field being 
the sum of the externally sweeped field W(i) and the 
molecular field from other particles that is determined 
self consistently. 23 This is a model of the nonlinear LZ 
effect that was applied to tunneling of the Bose-Einstein 
condensatei 2 ^* 2 ^ Again the problem can be reformulated 
in terms of a classical nonlinear LLE, see Ref.0, although 
the LZ effect is a quantum-mechanical phenomenon! The 
MFA solution shows that ferromagnetic interactions sup- 
press transitions since in this case the total effective field 
changes faster than W(t) and thus the effective value of 
the sweep-rate parameter e in Eq. (J5J becomes smaller, 
while antifcrromagnetic interactions enhance transitions. 
There are, however, models with couplings of different 
signs for which contributions to the molecular fields can- 
cel each other and the MFA wrongly predicts no effect of 
interactions. 

An important special case of Eq. JfjJ) is the so-called 
"spin-bag" model with the same coupling J for any two 



spinsA 2 ^ In the limit N -> oo with J = (N - 1) J = 
const the MFA for the spin-bag model becomes ex- 
act. This model can be exactly mapped onto the 
model with the spin S — N/2 and the Hamiltonian 
H = -W(t)S z - AS X - 2JS 2 . In contrast to the model 
with a general interaction that is difficult to solve numeri- 
cally, the spin-bag model can be numerically solved up to 
pretty big values of N that allows to check approaching 
to the mean-field limit N — > oo for the original model and 
to the classical limit for the equivalent spin-S" models A 
surprising result of Ref. that also should be valid for 
a general interaction is that MFA becomes exact in the 
linear order in small J, even without going to the limit 
N —* oo. Quantum corrections to the dynamics of the 
spin-S" model in the limit S 3> 1 (i.e., deviations from 
the mean-field results for the original spin-bag model) 
were systematically investigated in Ref. l2a 

For models with realistic finite-range interactions that 
are not small, the applicability of the MFA to the de- 
scription of the LZ effect is not justified. In particular, 
one can expect large deviations from the MFA for com- 
peting interactions such as the dipole-dipole interaction 
(DDI). Fortunately, in the fast-sweep limit £ < 1 one 
can construct a rigorous perturbative expansion of the 
staying probability P for Eq. © in powers of A 2 . The 
result can be written in the formSi 

P = 1 - e + e 2 /2 + e 2 I = e~ e + e 2 I (7) 

with Iq depending on the coupling J^-, resonance shifts 
Vi and the sweep rate v. The applicability of this expan- 
sion requires elo < 1. Eq. (JJJ shows whether interactions 
suppress transitions (Iq > 0) or enhance them (Iq < 0) 
that is not clear in the case of the DDL On the other 
hand, this rigorous result can be used to test the appli- 
cability of the MFA and other approximations that can 
be suggested to describe the many-body LZ effect. The 
results of experiments measuring P can be parametrized 
with the help of the effective splitting A calculated 
from Eq. |JSJ: 

A- = ^i, (8, 

as was done in Ref. 0. In the fast-sweep limit A eff should 
have the form following from Eq. @: 

Thus if Iq w const, as is the case for strong couplings, 27 
Eq. allows to determine A and Iq from the linear 
extrapolation to e — > 0, see Fig. [SJbelow. 

The aim of this article is to explain the derivation of 
Eq. J7J in a more detail, provide a comparison with the 
MFA result, and investigate the role of the DDI in crys- 
tals of molecular magnets of non-ellipsoidal shape where 
the magnetostatic field is inhomogeneous and thus the 
resonance fields for molecules in different parts of the 
crystal are different. The latter is needed to make a 
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comparison with the results of Ref. 11 on a crystal of 
a rectangular shape. 

The remainder of this paper is organized as follows. 
In Sec. [H]we construct the perturbation scheme for the 
many-body LZ effect at fast sweep and derive the general 
expression for Jo. In Sec. Illll we analyze different limiting 
forms of Iq that will be used below. In this section we also 
consider the role of static disorder, mainly due to nuclear 
spins, described by V. in Eq. © and perform averaging 
over Gaussian distribution of Vi. In Sec. II VI we consider 
the influence of the DDI on the Landau-Zener transitions 
in samples of the ellipsoidal shape, where the dipolar field 
is homogeneous, as well in samples of a general shape in 
the case of strong static disorder. We compare our results 
with the experimental data of Ref. [ll] in Sec. E] In Sec. 
IVII we recollect our main results and make some proposals 
for future research. 



II. PERTURBATION THEORY FOR FAST 
SWEEP 

We consider the transverse-field Ising model, Eq. JBJ, 
with the time-linear sweep W(t) — vt. The wave function 
of the system of N tunneling particles (TLSs) can be 
written as the expansion over the direct-product states 



*(*) = 



C 



mi , . ...niN 



(t)tp mi ®...®1p r , 



mi ,...,mjv — — 1,1 



1>-i 



1 Ih ^ 



I T). (io) 



The initial condition for is |C-i,...,-i(— oo)| = 1 

whereas all other coefficients are zero, i.e., the system 
starts in the state with all spins down. One can denote 
this state as \[[ . . . [) . As the time evolves, the state of 
the system becomes a superposition of all possible basis 
states in Eq. (|10fl . One can write the wave function in 
the form 

*(*) = C (t)|U...|)+^C i (t)<7 i+ m...|) 

i 

+|X> j (*)<7 i+C T J - + m... D+... (ii) 

ij 

including all possible numbers of flipped spins. Then the 
initial condition becomes |co(— oo)| = 1 and Cy(-oo) = 
0, etc. The normalization condition for this wave function 
is 



i = br + 



(12) 



The probability for a particle at the site i to stay in the 
initial state is given by 



\ c o\ 2 +^2\c 3 \ 2 + ^ 

3#i ' j^i,k^i 

i-M 2 -5>/ + ... 



(13) 



where we have used Eq. (|12|) to simplify the expression. 
The staying probability averaged over the system is 



p=-ypi = i--y\ci\ 2 --y 



(14) 



The Schrodinger equation for the coefficients in Eq. i|ll|) 
reads 

A 



ihc = x c - — } j 



A 



ihci = Ei(t)ci - —c - — ^ dj 

3 

ihdij = Eij(t)cij - ^ (cj + cj) - ^ y Ciji, (15) 



etc. Here Ei and Eij are the eigenvalues of the Hamil- 
tonian with A = and the ground-state energy E (t) 
subtracted 

Ei(t) = -W(t) + V h 

E t ,{t) = ^Wty + Vi + Vj-AJij, (16) 



where 



v; ; = K: + K (int) , v} int) = 



(17) 



Whereas Vi that enters Eq. © can be random, there is 

another contribution into resonance shifts, Vf nt) , that 
can gradually change across the sample for long-range 
interactions such as the DDI. 

For fast sweep there is little time for spin flipping and 
the system remains near the initial state: |co(t)| = 1 
while all other coefficients are small. Since spin flipping 
is caused by A, one can consider A as a formal small 
parameter and obtain the solution of Eq. \ib\ in the fast- 
sweep limit iteratively in powers of A. To this end, it is 
convenient to introduce slow amplitudes c(t) according 
to 

Cv (t) = ~c„{t)e- l *» {t \ <$>At) = ~ / dt'E v {t') (18) 

ft Jo 

with v = and rewrite Eq. 115(1 in the form 

dt c * ~ on 2^ c * e 



1 V- < .^x 



cue 



2h ^ ~ n 

3 



i[* 4 (t)-*y(t)] 



dt 2h 
i.g.. - i^g-e^-^W+^Wl + (i — i) 

dt 13 " 2h Cl6 



(19) 



In the last equation we have dropped the term with Cy; 
since we are going to calculate P up to A 4 and it can be 
shown that at this order Ciji is irrelevant. Iterating Eqs. 
(EH yields 

co = 1+(|) 2 4 2) + ... 
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g. - !^c (1) + ( —) c (3) 



!AV (2) 
2h ij 



(20) 



where we have written down only relevant terms that are 
given by 



J — CO 
t 



— oo 

t 



1 

ft 



£ / di'cg^iOe^*^*')-**^*')]. (21) 

J — oo 



Now the expansion of Pi of Eq. (|13[) up to A 4 has the 
form 



Pi(t) = 1- 



A 

2ft 

E 



c, (1) 



_(2) 
r- ■ 



2 Re 



(c«*cf>) 



.(22) 



It is convenient to split Pi into the part corresponding to 
noninteracting particles and the part depending on the 
interaction, the result for the former being already known 
from Eq. (JSJ, and to introduce e. Thus one obtains at 

t — ► oo 

P l ^l-e + e 2 /2 + e 2 J2 h - e~ e + e 2 ^ J« (23) 



with 



V27rft/ 



2Re<!cP 



(oo) 



<V)W {t)e -Mf(t) _ gg) (t)e -<*«(t)] | (24) 



Here the quantities with the superscript (ni) correspond 
to the noninteracting system, Jij — 0. 

There is a seeming paradox in the derivation of Eqs. 
(|23|l . I|24|) . For a macroscopic system, TV 3> 1, even at 
fast sweep, e <C 1, a macroscopic part of the particles eN 
is flipping out of the starting state with all spins down. 
However the wave function of Eq. lllfl with only one 
or two spins flipped was used in the calculation above, 
whereas the coefficients Cijk etc. have been dropped. Our 
calculations is nevertheless correct, because extending of 
the expansion of Eq. I|2U|I to include more terms would 



only lead to contributions of orders higher than A 4 or, 
correspondingly, higher than e 2 that we are not inter- 
ested in. The seeming paradox has nothing to do with 
the interaction and it emerges already for a system of 
N identical noninteracting TLSs. In this case one can 
just solve a one-particle problem exactly and expand P 
of Eq. JSJ up to e 2 or, alternatively, one can solve the 
one-particle problem perturbatively using lines 1 and 2 
of Eq. (HI. The result is P = 1 - e + e 2 /2 + .... On the 
other hand, one can do the perturbative expansion for 
the whole system of N particles up to e 2 or higher that 
evidently leads to the same result. The key observation is 
that at any order in s (both with and without the inter- 
action) the terms that are diverging in the limit N — *■ oo 
cancel each other, as it should be. In particular, one can 
see from Eqs. ()23|l . I|24|) that there is no proplem in the 
thermodynamic limit since Iij decays with increasing the 
distance between i and j if the interaction does so. 
In connection to the above it should be mentioned that 
the truncated Schrodinger equation, Eq. 1|19[) cannot be 
directly solved numerically. The problem is that the so- 
lution of this equation contains all powers of A or e, and 
starting from A 6 or e 3 not all terms of the same order 
are taken into account. This leads to noncancellation of 
terms diverging in the limit N — > oo (i.e., spurious terms 
of order e 3 N etc.) and thus to an unphysical behavior of 
the whole solution. 

Calculation of the double and triple time integrals in 
Eq. H24fl requires essential efforts. At least for the aver- 
age staying propability of Eq. I|14|) where only the sym- 
metrized quantity 



r(sym) 



1 



{Iij + Iji) 



-J — *% 3 2 
is needed and from Eq. Ijl4(l one obtains Eq. (JJJ) with 



(25) 



if sym ' can be expressed as 



r( s Y m ) 

«3 



(26) 



I. 



(sym) 



Ai 



(27T%/3 y ) By + sin (27r%/3 i3 -) 



ij ■ 

Here Aij , Bij , are combinations of Fresnel integrals 
C(x) and S(x) 



1 i 

2 ~ 4 



C(7«) 



s ha) 



1 


"1 




2 i 




4 


2 


-cha) 


4 


\ - s hjd 



(28) 



Bi 



cha) 

s ha) 



■Cha) 
s ha) 



(29) 



2-cha) 



2 - s ha) 
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FIG. 3: /f! ym ' vs its independent arguments ft ■ and Si 



FIG. 4: J^ ym) as a function of the resonance shift Sij for 



ft, = ±5. 



2 - 5 ( 7 »j) 



2- C (7ii) 



(30) 



and the dimensionless variable 7^ is defiried by 



7ij = + 
- = Yi 

a i — 1— — — ! 

vzirnv 



S i: j n, OCj 



Pi 



4Ji. 



4J f 



V2irhv ttA 



/e. (31) 



Here Vi contains interaction according to its definition in 
Eq. I|17fl . Note that independent arguments in Eqs. J57J- 
ffST7f> are reduced interaction (3^ and resonance shifts 6^ , 
thus we also will be using explicit notations of the kind 



(32) 



for any of the functions , Ay, -By, and Cy. S 

larly to Eq. (|17|l . one can write 



<5y = <5y 



c(int) 



(33) 



Eqs. and l|26|) - (|31[) is our main result that is valid 
for arbitrary interactions Jy and resonance shifts Vi. 
Note that it has a pair structure and thus it can be ver- 
ified against the direct numerical solution for the model 
of two coupled particles. Analytical form makes its ap- 
plication practically possible: Triple time integrals of Eq. 
(|24|l cannot be computed numerically with a reasonable 
precision within a reasonable time. 



III. ANALYSIS OF THE SOLUTION 



A. General properties and limiting cases 



A three-dimensional plot of -/y ym ^ vs its independent 

arguments /3y and Jy is given in Fig. [3J One can see 

that I^ ym ^ has a plateau ly = 1/2 for strong ferro- 
magnetic interactions, (3.^ 3> 1. This means, according to 
Eqs. J7J) and l|26[l ferromagnetic interactions suppress LZ 



transitions. For strong antiferromagnetic interactions, 
—fiij ^ 1, the value of I^™ 1 ^ oscillates between the upper 
bound 1/2 and the lower bound around — 3/2(see Fig.0J, 
i.e., on average antiferromagnetic interactions enhance 



transitions. For large resonance shifts, 



'j 



Z \0ii\ , the 



value of iy decays to zero. This "causality" is phys- 
ically expected since two particles having resonances at 
very different values of the sweep field W(t) do not affect 
LZ transitions of each other. 



It can be shown that I, 



(sym) 



satisfies the sum rule 



dSI {sym \S,/3)=/3 



(34) 



(sym) 



that drastically simplifies the analytical results in the 
case of strong static disorder. 

In the homogeneous case Jy — one has /, 

Jy = 1(0, p tj ) = F(f3 tJ ), where 

F((3) = C ((3) [1-C (/?)] +S{0)[1-S ((3)} 

The limiting forms of F([3) are 



(35) 



-§-^cos(f/3 2 

F{J3) - { - (3 2 , 

1 _ 1 

2 7^9)*' 



5), 

\P\ « 1 
/?> 1. 



(36) 



For the weak interaction, |/3y | <C 1, Eq. (JJJ then yields 
at the leading order 



P = 1 — e 



4Jo 5/2 

7rA ' 



(37) 



a generalization of Eq. (26) of Ref. for the arbitrary 
form of Jij. Note that Eq. (|37() is essentially a MFA re- 
sult (see a separate consideration of the MFA for our 
model in the Appendix) as it only depends on the zero 
Fourier component Jo of the coupling Jy. In contrast 
to thermodynamic systems, here the applicability of the 
MFA is controlled by the strength of the interaction in 
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addition to its radius. The correction to the mean-field 
result of Eq. I|37[l is described by the term — (3 2 in the 
central line of Eq. (|36[) . For the nearest-neighbor inter- 
action with z nearest neighbors the relative correction 
to the last term of Eq. is - [4J /(ttA)] e 1/2 /z that 
becomes small both for large z and for fast sweep. 

The saturation for strong ferro- and antiferromagnetic 
interactions in Eq. (|36[> corresponds to the case of well- 
separated resonances studied in Sec. Ill of Ref. |!J At 
fast sweep e <C 1, LZ transitions happen in the range 
W ~ Wlz around the level crossing, where 2 ^ 



W LZ = V2nhv = nA/y/e. 



(38) 



That is, in Eq. (jHU (3 i3 = AJij/Wlz- The resonances 
are well separated for |Jy| 3> Wlz, i.e., for |/3jJ ^> 1. 
This limit cannot be described by the MFA. The latter 
becomes valid in the limit of nonseparated resonances, 

Let us proceed to the inhomogeneous case, Sij ^ 0. 
For 1/3^1 < 1 Eq. yields 



(39) 



with 

m 



in (|«5 2 ) + cos (|<5 2 ) + nS [S (S) - C (5)} 



V2 



ir8 2 /2, 

(¥ 2 + 



f) 



S 2 « 1 
S 2 » 1. 



S(5) satisfies the sum rule 

/oo 
d<5S(5) = 1 
-OO 



(40) 



(41) 



that is a particular case of the more general Eq. 134JI . We 
will see in the Appendix that Eq. (|39|) also follows from 
the MFA. t 

max(|/9y|, 1) > 1, Eq. yields a small 



For 
value 



r(sym) 



V2 



■ cos 



(1 2 



+ 



Pi 



6 >3 



Pi 



(42) 



as explained above. Note that here r y ji = —6^ + (3^. For 
5ij ^> max(|/3y |,l) this reduces to Eq. (|3~9"j) with the 

second limiting form of Eq. I|4l)[l . That is, the applicabil- 
ity range of Eq. (|3~9l is larger than just -| <C 1. 

In the case • | - 



^S> 1 one obtains 



f (sym) ^ j -1/2 - cos \2iT~8ijP i ^ , (3 t 



< 



1/2, 



0« > 0. 



(43) 



For long-range interactions such as the DDI each TLS can 
strongly interact with many other TLSs with different 
strenghts 3> 1, and the resonance shifts between 

different particles can be strong and different, \5ij \ 1. 



In this case the cosine term in Eq. (1431) averages out in 
Eq. H§. 

Neglecting the small value given by Eq. I|42(l and re- 
placing the cosine term by zero in Eq. (|43|l one can com- 
bine the expression for I^™ 1 ^ in the case of both large 
arguments and 5^: 



I 



(sym) 



sign(/y (l/y -|%|). (44) 



Here 9{x) is the step function. Note that this form sat- 
isfies the sum rule, Eq. (|3~ljl . 

The terms with cos (2n~& \j ft and sin (^irSijP^ in 

Eq. (|27|l represent the effect of the quantum-mechanical 
phase in the many-body LZ effect. The occurence of 
these terms is due to the possibility to come to a given 
final state along different ways. The latter can be seen 
in Fig. [21 but is absent for the usual LZ effect, Fig. ^ 
The quantum-mechanical amplitudes corresponding to 
the different ways add up with their phases in the ex- 
pression for the staying probability P that causes its os- 
cillations. The minimal model that exhibits this effect 
is a dimer of two antifcrromagnetically coupled two-level 
systems with shifted resonances^ Note that the effect 
of the quantum-mechanical phase is totally absent in the 
mean-field approximation. 



B. Effect of static disorder of resonance positions 

As can be seen from Eq. l|17f) . the shifts of the reso- 
nance positions Vi that enter the final formula via Eq. 
(|31|l are the sums of the two terms: The original shifts 
Vi of Eq. © and the contribution of the interaction Jy . 
The former can arise due to different types of static dis- 
order, including that induced by nuclear spinsia (see 
below). The latter are constant and thus irrelevant 
^(mt) _ y^(mt) ^ ^ w ^j^ n ^- ne body of the sample for 
short-range interactions, and in this case they only affect 
the particles on the surface. For long-range interactions 



such as the DDI, V 



(int) TA (int) 



V in general smoothly varies 



across the sample, depending on the sample shape. 

The effect of static disorder can be accounted for 
by averaging Eq. (|27|) over stochastic values of on — 
[Vi j (7rA)] y/e with a normalized Gaussian distribution 

1/2 

p a {a) = (2na 2 ) exp [—a 2 / (2a 2 )] and quadratic av- 



erage 



<«?> 



a . The distribution of Si 



a, 



a,- is 



then given by the same function with a 2 
ing over the static disorder 



2(7 . Averag- 



7 (sym) (a,5 (int) ,/3) 



dS 



: exp 



) V Ana* 
xI ( - Bym) (S + .) 



'4a 2 

(int 



[3) (45) 



can only be done numerically in the general case. There 
are, however, particular cases in which it can be done 
analytically. 

In the case of strong disorder a ^> max(l, |/3|), the 
integral in Eq. I|45[) converges at S <C er, so that one can 
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neglect the exponential, shift the integration variable and 
use Eq. J32JI that yields 



I 



(syn 



cr ' 



a » max(l, 1). 



(46) 



Note that in this limit I^ ynl ^ is independent of Sfj . For 
— and a 1 one can use the simplified form of 
/| sym ' ) given by Eq. l|33Jl to obtain 



f(sym) 



■ erf 



2(7 



CT > 1. 



(47) 



Indeed, for ct ~ \ fi i A >• 1 Eq. (|l4*jl becomes correct while 
reducing leads to Eq. (|47)|l . 

In the case \(3\ <C 1 one can use Eq. {59). where the av- 
eraged value of S can be found analytically for arbitrary 
cr in the case J*- 111 *) = 0. One obtains 



ij ym) - ^s(ct), 



where S(cr) is given by 

H(ct) = /(2^ct 2 ), /(a) 
X + 1 



r5 (ir 



0. 



(48) 



i/2(l + x 
x - 1 



vVi + a^ + i VVi + x 2 



.(49) 



Here /(x) monotonically decreases and has limiting forms 



/(*) = 

The latter limiting form, 
H(ct) 



1 - x/2, x < 1 
l/V^x, x>l. 



20Fcr' 



(50) 



(51) 



follows from Eq. (|4Tj) and it is valid for <5^ lnt ^ 7- as well. 

Let us now consider the case of weak disorder a <C 1, 
setting (5' mt ) = 0. In this case the main effect arises for 
large negative /3 [see Eq. (|43f) ]. where even small disorder 
can average out fast oscillations of I^™^ ■ Combining Eq. 

J2J with Eq. g5J and using (cos (2tt<5/3)) = e -( 2 ^) 2 one 
obtains 



r(sym) 



F^ + ei-P^^l-e-i 2 ^) 2 ). (52) 



It can be seen that the limiting value lim^^_oo F((3) -= 
—3/2 in Eq. I|36[) is unstable with respect to 
small static disorder: For any a > one obtains 
lim^-oo I( sym V,0,/3) = - V 2 - Th e reason for that is 
the effect of the quantum-mechanical phase that leads to 
the fast oscillating terms in Eq. (|27() . 

/^ sy averaged over static disorder is shown as a func- 
tiion of (3 for different a in Fig.[S] One can see that a 7- 

-co 



changes the asymptotic behavior of I, t 



(sym) 



at /? 



0.5 
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FIG. 5: /^ ym ' as a function of fi for different strengths of 
static disorder a. 
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15 



and that i^ sym ' ) becomes odd in (3 for a 1 (in fact, 
already for cr > 1), as described by Eq. I|47|l . 

In molecular magnets distribution of resonance posi- 
tions for electronic spins V.- L is mainly induced by the 
coupling to nuclear spins. This can be both the contact 
hyperfine interaction with the nuclear spins of the mag- 
netic atoms, as in Fes containing the isotope 57 Fe and 
in Mni2 , and the dipole-dipole interaction with magnetic 
moments of nuclear spins of non-magnetic atoms, mainly 
protons. The latter is small as the nuclear magneton but 
the number of nonmagnetic atoms interacting with the 
magnetic atoms is of order 10 2 , so that the total effect 
is quite substantial. Measurements of Ref. for the 
resonance between the ground states \±S) (S = 10) in 
Fes with the standard Fe isotope (no nuclear spins on 
Fe atoms) yield a Gaussian line shape with the width 
cth — 0.8 mT that is in agreement with the theoretical 
evaluation gh = 0.5 mT«i4 For the 57 Fe isotope the line 
width is only about two times larger, 1.2 mT and 1.1 mT, 
respectively^ The dimensionless dispersion a introduced 
earlier in this section is defined similarly to Eq. (|31[) : 

(J = /== = ^oVe, (53) 

where cry = 2Sgfi B an is the dispersion in the energy 
units. With <th = 0.8 mT one obtains cry /fee = 21 mK. 
With A/fcs ~ 10~ 7 K for the ground-state resonance one 
obtains a huge value cto — 0.68 x 10 5 that makes cr very 
large even for a moderately fast sweep such as e ~ 10 ~ 2 . 
With the help of Eq. (|31|l one can rewrite Eq. I|46[) in the 
natural form as 



(sym) 



2 Jij 
sf^ 'cry 



\JiA < cry 



(54) 



that is independent of A. One can show that this result 
follows from the MFA as well, as all our limiting cases 



where /, 



(sym) 



Jij. If the disorder is so strong that 



| | -C cry is satisfied for all distances, including the 
nearest neighbors, then the many-body LZ effect can be 
described by the MFA. 
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It remains only to justify that nuclear spins can be 
considered as static disorder. The appropriate condition 
is 



(55) 



where ujlz the characteristic frequency of the LZ tran- 
sition and r, mc is the nuclear relaxation rate. At fast 
sweep, £<1, LZ transition occurs in the range of sweep 
fields defined by Eq. I|38|) . This yields the Landau-Zener 
time thz = Wlz/v = \J2+h/v and the LZ frequency 



ulz = 



tLZ 



V 

2Trh 



A 1 
2h~Je' 



(56) 



For A/k B ^ 1CT 7 K one obtains uj lz = 0.65 x 10 4 /Ve 
s _1 that grows with the sweep rate. On the other hand, 
recent measurements 3 ^ on Mni2 yield F nuc in the range 
between 10 -1 and 10 -2 s _1 for T < 1 K. Thus nuclear 
spins are really slow and they can be considered as static 
disorder in a wide range of sweep rates. The situation is 
unlikely to be strongly different in Fes- 



IV. DDI IN CRYSTALS OF ELLIPSOIDAL 
SHAPE OR WITH STRONG DISORDER 



(~ 600 A for Fes) that strongly interact with a given 
spin, \(3, tj \ > 1. 

As can be seen from Eq. (|46l) . in the case of strong 



static disorder, 

reduced if | 3^ | 
length is 



er 3> 1 the value of I, 



(sym) 



is strongly 



< 



a. The corresponding characteristic 



(61) 



Within this distance the interaction is still strong, 
1. The estimation for the number of spins 



rOy m ) 

ij 

stronly interacting with a given spin is 
- 3 £ 4E D 



v a 



(62) 



For Fes one obtains N a ~ 24 that is large but much 
smaller than the values of order N% ~ 10 5 in the absence 
of disorder. For the superstrong disorder, 



a > £ or cry > E 



D, 



(63) 



< 



there is no range where the interaction is strong, r a 
v^ 3 . For ay ^> Ed, as commented after Eq. (|54H . the 
mean-field description of the LZ effect with interactions 
becomes valid. 



Let us now turn to the DDI between tunneling spins 
±5 of magnetic molecules aligned along the easy axis 
pointing in some direction m, |m| = 1: 

3 cos 2 6ij - 1 . . 

Jij=E D <Pij, 4>ij=vo 3 , (57) 

where 

E D = (g^ B S) 2 /v (58) 

is the dipolar energy, vq is the unit-cell volume, r\j is 
the distance between the sites i and j, and cos 6 \j = 
(ry m) /rij. The two mostly well known molecular mag- 
nets are Mni2 and Fes, both having the effective spin 
S = 10. Mni2 crystallizes in a tetragonal lattice with 
parameters a = b = 17.319 A, c = 12.388 A (c is the 
easy axis) and vq = abc = 3716 A 3 . Fes has a triclinic 
lattice with a = 10.52 A (a is the easy axis), b — 14.05 
A, c = 15.00 A, a = 89.9°, (5 = 109.6°, 7 = 109.3° and 
vq = abc sin a sin (3 sin 7 = 1971 A 3 (see, e.g., Ref . l3l|) . 
Fes is more convenient as a model system for us as (i) 
the standard Fe isotope does not have a nuclear spin that 
is ignored in our theory and (ii) the splitting A in Fes 
has a well-defined origin and it can be estimated theoret- 
ically. One can write f3 { j of Eq. I|31|l in the form 

/3fc- = e^, ^V 1/2 , (59) 

ForFe 8 E D /k B = 126.4 mK and A/k B ~ 10" 7 K,sothat 
£q ~ 1.6 x 10 6 . Thus for not too fast sweep, e ~ 10~ 2 , one 
has £ ~ 10 5 . This is also an estimation for the number of 
spins N% ~ ^ within the distance 

r c = (vodf 3 (60) 



A. Samples of ellipsoidal shape 

Consider a macroscopically large specimen of ellip- 
soidal shape. In this case the magnetostatic field inside 
the homogeneously magnetized sample (the system re- 
mains in the vicinity of this state in the LZ effect at 
fast sweep) is homogeneous in the bulk of the sample, 
£(mt) _ q Tj-^g one only has to make averaging over 
static disorder using Eq. |@SJ|. According to Eq. l|4"5|) 
j(sym) j oeg na |. (jjygj-gg f or \a I ^ 00. For $ 3> 1 and 

<y <C C (i- e -7 A'cr ^> 1) one can replace the sum in Eq. I|26() 
by an integral converging at r = ~ min(r c ,r (T ) that 
makes the result independent of the lattice structure: 

Io= f— I (sym) (a,/3(r)), (64) 

J Vq 

where 

/3(r) = ^(r), 0(r) = wq 3008 ^" 1 . (65) 

To simplify the integration in Eq. (|64|l , it is convenient 
to change the variables and integrate over the direction 
r/r and the value of /3 instead of integrating over the 
direction r/r and the distance r. At large distances where 
(3 <C 1, one can use Eq. (|48|l that makes account of the 
static disorder. Since 7 < - sym ^(cr, /3(r)) oc /3(r) behaves as 
the DDI, the result of the integration depends on the 
sample shape. Thus one has to be cautious with changing 
integration variables. One of possible ways to tackle this 
problem is to introduce a subtraction function, say 

F subtr (a,/?(r)) = f (r 2 } . . S(a) (66) 
1 + 8 (r) 
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that has a sufficiently simple form, does not diverge at 
small distances, (3 — > oo, and has the same behavior as 
li/ at large distances, |? « 1. (Other types of sub- 
traction functions differing by the cutoff at (3 — ► oo, for 
instance, the function having 1 + |/?(r)| in the denomina- 
tor, yield the same results.) With the help of Eq. jffijjl 
one can write the integral in the form 



h = 



d 3 r 
vo 

+ 



J^ m )(a,/3(r))-^ subtr (a,/3(r)) 

/ — F subtr (a,/3(r)). (67) 
J v 



The first integral converges fast at large distances, thus 
it does not depend on the shape and can be rearranged 
by changing variables as said above. After doing that, 
the contribution of -Fsubtr into this integral disappears be- 
cause of the antisymmetry F su btr(-/3, c) = -F su ur((3, <?)■ 
The second integral can be calculated analytically for the 
actual sample shape using the results of the magnetostat- 
ics, without changing variables. The final result has the 
form 



Iq == G£, G — Gj + G 



Here 



subtr 



Gi = %v r d/^t^ 



9VS 



(68) 



(69) 



where I^ sym \a,(3) is given by Eq. ll4*5|l . The subtraction 
contribution reads 

G subtr =\k+ fl/3 - n^) 4tt1 H(ct), (70) 



where the constant K is given by 

K^-^(l-±\n ^±l) = -0.66924, (71) 
9 \ V3 V3- 1 J 

and 

n (m) ig 

the demagnetization coefficient depending 
on the sample shape and on the direction of the easy 
axis (i.e., the magnetization) m. If m coinsides with the 
z axis, the symmetry axis of the ellipsoids of revolution, 
then n^™- 1 = = 1/3, 0, and 1 for a sphere, needle and 
disc, respectively. In Eq. I|70|l the term with K is the 
result of the actual calculation for the sphere, whereas 
the remainder is known from the magnetostatics. 

For a <C 1 we use Eq. (|52[) to calculate G/. One obtains 
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FIG. 6: G of Eq. I|68|l for the sphere vs width of distribution 
of individual resonances a. Dashed lines on the left and right 
are asymptotes of Eqs. I|72|l and l|73|l . respectively. 



the leading term in G is G sub tr- From Eqs. ij7U|) and l(5T7|) 

one obtains 



G^ 



K + (1/3- n( m )) 4tt 



2^a 

Using Eqs. lO, EH, and (O, one obtains 



In ^ 



K + (1/3- n (m) ) 4tt 



E D 



(73) 



(74) 



in the case of strong static disorder. For the ground- 
state resonance in Fes the values of I are 23.4, —4.44, 
—60.0 for the needle, sphere, and disc, respectively. Io 
changes its sign for the critical value of the demagnetizing 
coefficient defined by 



(75) 



One has Iq > for m m ) < n 



(m) 



and Iq < for n( m) > 



The results of the numerical evaluation of G for differ- 
ent values of a are shown in Fig. Whereas G < for 
the sphere and disc, DDI acting predominantly antifer- 
romagnetically and enhancing LZ transitions, the result 
for the needle in Eq. (|68|l becomes positive already for 
a > 0.1. 



G = — 5 734 + 16 (7r/3) 5 / 2 a + f 1/3 - n (m ^ A-k (72) B- Superstrong disorder - the mean-held limit 



Here the linear-cr contribution stems from the second 
term of Eq. 152|) that changes the asymptotic behavior of 
/' sym '(<7,/3) at (3 — » — oo. The large numerical factor in 
front of this contribution makes G very sensitive to a. 

For a 3> 1 the integrand of in Eq. (|69|l is of order 1/a 
and, in addition, it is nearly odd in (3. Thus the result of 
the integration is smaller than 1/a. In the limit a 3> 1 



In the MFA, I^ ym) is given by Eq. that yields 
Eq. H48fl after averaging over static disorder. We have 
seen above that the MFA becomes valid for the descrip- 
tion of the LZ effect only in the case of weak interac- 
tions \Jij \ <C A/^/eor in the case of superstrong disorder 
ay 3> max(| Jij |) if the interaction is strong. In molecular 
magnets the applicability condition for the MFA is not 
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fulfilled. Yet we consider it for the sake of completeness. 
For the DDI one obtains 



r(MFA) _ 1 C n 2 E r 



-D„ 



(76) 



V— ~ mm i— 

where D mm = m • D m 

D mm = £ ^ = D£ m) + (l/3 - n< m >) 4vr (77) 

j 

and D m is the dimensionless quantity proportional to the 
dipolar field 



9VbS 



vo 



(78) 



created on the lattice site i by all other spins, j =/= i, 
pointed in the m direction. One can see that D m m re- 
places the universal constant K in Eq. (|74|l . 

For ellipsoids of revolution with magnetization m di- 
rected along the symmetry axis, D mm is independent of 
the lattice site i in the main part of the sample, except for 
the vicinity of the boundaries. Note that the macroscopic 
field theory (magnetostatics) is insufficient to obtain Eq. 
(|77|l . To this end, one can introduce a macroscopic sphere 
around the site i. The field from the spins at sites j in- 
side this sphere can be calculated by a direct summation 
over the lattice and it yields the first term in Eq. (|77|l . 
The field from the spins outside this sphere can be calcu- 
lated macroscopically and it results into the second term 
in Eq. J77J). For a simple cubic lattice fi c rc ' = by 
symmetry and the result for D zz becomes purely macro- 
scopic. For tetragonal lattices -Dmm Clc ' ) >0 if a = b > c 
and fl£ clc ' < if a = b < c. Direct numerical cal- 
culation yields I>Sm ere) = 5.139 for Mni 2 and 4.072 for 
Fes- Note that Eq = — (1/2)D zz Ed is the dipolar energy 
per site for the ferromagnetic spin alignment. Our result 
Eq = — 4.131.Ed for the needle-shaped Fe« is in qualita- 
tive accord with E = -4.10E D of Ref.Efl. 

One can see that unlike the rigorous value of Iq in 

(MFA) 



the limit £ ^> 1, the value of Iq depends on the 
lattice structure via 

D (Sphere) Eq ^ ^ 

K flfSi here ' yields ni m) = 0.742 (i.e., 4 MFA) > for 
n ( m ) < 0.742), whereas for Fes one obtains ni" 1 ' = 0.657. 
In both cases it essentially differs from the prediction of 
our rigorous theory, Eq. 1751) . 

C. Samples of general shape with strong disorder 

The strong-disorder result of Eq. (|74|) is valid for the 
samples of nonellipsoidal shape as well. Indeed, the 
derivation of Eq. I|74|l is based on the sum rule, Eq. 
(|34|l . that makes the integrand in Eq. i|69|) linear in (3. 
This alone would be, however, insufficient, as the limits 
a — ► oo and (3 — ► oo (i.e., r — > 0) are not interchangeable. 
Eq. H46f) diverges at r — + and it needs an appropriate 

regularization, whereas the true ly does not diverge 
at r — > 0. The regularization we have used above consists 



in choosing j3 as the integration variable. Then approxi- 
mate oddness of ly m ' in (3 for large a strongly reduces 
the integral over f3 and makes it of order 1 /a 3 that can 
be neglected. This oddness is also preserved in the limit 
(3 — > oo for large but fixed a, as follows from Eq. 144fl . 
The result of the regularization at r — > is the nontrivial 
constant K in Eq. (|74|l . One can see from Eq. (|44|l that 

J^ sym ' ) remains approximately odd in (3 for large (3 and 
a even in the case 5^™^ ^ for nonellipsoidal shapes, 
especially as 8^ nt ^ — > for r — > 0. Thus the situation 
in the strong-disorder case a 3> 1 is the same for ellip- 
soidal and nonellipsoidal sample shapes. In both cases 
one switches to the integration over (3 and obtains zero 
because the integrand is odd. Recalling that integration 
over j3 is incompatible with taking into account the sam- 
ple shape, one can correct the situation by introducing, 
e.g., the subtraction function of Eq. (|66|l . 

The general result thus can be written in the form 



d 3 r 

Io- I — 



d 3 r' 5(<r)/3(r' - r) 
~l + /3 2 (r'-r)' 



(79) 



that contains the additional averaging over the sample, 
S(cr) being given by Eq. 1)5 l[l. For the spherical shape 
the result of the integration is known from Eq. (|74)l . Iq = 
£S(a)K. The difference between the general shape and 
the sphere arises from the integration at large distances, 
where (3 1 in the denominator of Eq. I|79|) can be neglected. 
Thus with the use of Eq. I|65|l one can write 



K 



F(r) 



(80) 



where the difference F(r) is given by 



F(r) = f flL^r'-r)- f — 4>(r> - r) 

^Sample v ./Sphere v 

= M')-^. (si) 

Here T> mm (r) is the m-component of the reduced macro- 
scopic internal field created by the sample magnetized in 
the m direction, 



mm V 1 J i 



(82) 



and £> mm 

yields 



(Sphere) 



(2/3)4?r. Inserting Eq. JST) into Eq. 



K + (l/3-n (m A 4tt 



E D 



(83) 



[c.f. Eq. lOJ], where the average demagnetization coef- 
ficient is given by 



1 f d 3 r^ , , 



(84) 



For the needle, sphere, and disc with m along the symme- 
try axis one has 2? mm = 47r, (2/3)47r, and 0, respectively, 
so that Eq. I|84() recovers the known results for these ge- 
ometries. 
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FIG. 7: Dimensionless macroscopic internal field X>mm of Eq. 
1891 in the sample of the rectangular shape with L x — L y — 1 
and L z = 3, uniformly magnetized along the z axis, plotted 
vs x and z for y — 0. 



For the general shape the integrals in Eq. (|81|l can be 
calculated directly with the use of 



(r) = m-VQ 



3n r (m • n r ) — e z 



m-^orot 



[m x n r ] 



where n r = r/r. The integral formula 



dVrotF 



dSxF 



(85) 



(86) 



for an arbitrary vector function F(r) then reduces the 
integration to the surface. Alternatively one can use the 
the Bio-Savard formula 



H(r) = - I d 



,j x (r-r') 



(87) 



with j =crotM and M — (g/j, B S/vo) m, |m| = 1. As the 
molecular currents j are nonzero on the surface only, one 
obtains 

H m (r) = ZtfS [ [dS' x m] x J^JL. (88) 
vo Js |r' - r| 

For the samples of the box shape with sides 2L X , 2L y 
and 2L Z with generally directed m the integration in Eq. 
(|88|l yields the sum of 24 arctan terms for the symmet- 
ric part of 2? mm (r) that makes the contribution into the 
volume average in Eq. (|54*}l 



arctan 



,-1 



{L x + T] x x) (L v + r) y y)(L z + r] z z) 

\J {L x + V x xf + (L y + T] y y) 2 + (L z + r) z zf 
+ (xyz yzx) + (yzx zxy) . (89) 

If m is directed along one of the symmetry axes of the 
box, the antisymmetric terms in T> mm (r) disappear and 



this formula yields T> mm (r). The result above is illus- 
trated in Fig. for m z — 1. Adopting Eq. lj59"|) in Eq. 
(|SH> yields = 1/3 for the cube and ~ 0.102 for 
the box with proportions L x : L y : L z = 1 : 5 : 8 that is 
close to the shape of the crystal used in Ref. 0. 



V. COMPARISON WITH EXPERIMENT 

In the sweeping experiments 11 studying the ±S tran- 
sitions in Fes the standard LZ effect can be seen down 
to dB/dt ~ 0.01 T/s. Using v = 2Sg^ B dB/dt and 
e = TrA 2 /(2hv) with A/k B ~ 10 -7 K one obtains the 
range e < 10 -2 for the standard LZ effect. In the re- 
gion of slower sweep dB/dt < 0.01 T/s (i.e., e > 10" 2 ) 
the effective splitting A cS calculated from Eq. JSJ) goes 
down. This suggests that here the Landau-Zener effect 
is strongly modified by interactions that suppress transi- 
tions. It is very instructive to replot experimental data^ 
as P vs e, see Fig. [Ht- The high plateau of P(e) at large 
e is similar to that for the spin-bag model with the fer- 
romagnetic coupling and large N, see, e.g., Fig. 3 of Ref. 
9. One can estimate lo if one replots the experimental 
data for A cff vs e using Eq. see Fig. [5] One can see 
that A cff is apparently linear at small e. The fit ignoring 
the downward bump near the origin yields /q XP ~ 11. 

The shape of the crystal used in Ref. [llj was an elon- 
gated platelet (l z = 80 /im, l y = 50 fim, l x — 10 fiwM>). 
The crystallographic easy axis a slightly deviates from 
the direction of the longest axis of the crystal, l z . It is 
rotated by a ~ 9° in the yz plane and then rotated by 
(3 = 7° away from the yz planed Thus for this shape 



m z = cos a cos (3 ~ 0.9803 
rriy = sin a cos /3 ~ 0.1553 
m x = sin/3 ~ 0.1219. 



(90) 



Neglecting this small tilt (i.e., setting m z = 1, m x = 
m y = 0) and using Eqs. (|84|) and (|89ll . one obtains 
7l (m) = n (z) _ o.l02. Then Eq. ® with E D /k B = 126 A 
mK and <Jyjk B = 21 mK yields Iq ~ 14.8 that is in a 
good accord with 7q XP ~ 1 1 . Taking into account the tilt 
of the easy axis improves the agreement with the exper- 
iment: n< m ) ~ 0.113 and J ~ 13.8. 



VI. DISCUSSION 

We have analytically investigated the asymptotic stay- 
ing probability P in the many-body Landau-Zener effect 
at fast sweep, e< 1. The coefficient Iq in the expansion 
of Eq. (J7J) accounting for the interaction has been cal- 
culated rigorously for arbitrary interactions J,j and the 
resonance shifts Vi in the transverse-Ising Hamiltonian 
of Eq. ijlj). We have shown that ferromagnetic interac- 
tions (J > 0) increase P (i.e., suppress LZ transitions) 
whereas antiferromagnetic interactions ( J < 0) act in the 
opposite direction, for the fast sweep. 

The resonance shifts Vi [see Eq. I|17|l ] have been shown 
to reduce the value of Jq since different particles under- 
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field reduce the value of Iq. The tunneling particles effec- 
tively decouple at the distances satisfying • VD ZZ ~ 
(f)^ . Estimating the gradient as VD ZZ ~ L, where L is the 
linear size of the specimen, one obtains the characteris- 
tic length r; n h = (Lvq) 1 ^ 4 ■ In the strong-gradient case 



< 



(V/vq) 1 / 4, , with V ~ i 3 . This yields still large values 
of Iq for macroscopic samples. Jo even increases if the 
sample is elongated or flat because in this case the field 
becomes more uniform than in the cube. One obtains 
In ~ — 10 4 for the crystal used in LZ experiments of Rcf. 
fill if one neglects the static disorder. 

Taking into account interactions with nuclear spins 
considered as frozen-in disorder drastically reduces Iq for 
Fes and yields values that agree with the values extracted 
from the measurements of Ref. ^] both in sign and mag- 
nitude. 

It would be important to perform LZ experiments on 
crystals with a more flat shape for which the demagnetiz- 
ing coefficient > 0.28 and Eq. (3J yields I < 0. In 
this case DDI enhances LZ transitions, and the value of 
A e g increases with e starting from its initial value A, c.f. 
Fig. |SJ An interesting question is whether this tendency 
holds for larger e (slower sweep rates) as well. Our theory 
is applicable only for small e and it cannot answer this 
question. It is possible that A e g has a maximum at some 
e and then it falls below A, as indicated by our results 
for some simplified models of the many-body LZ effect. 



one can use Eq. (|44|l that yields Iq ~ rf nh /vQ 



o.oo 



0.05 



0.10 



0.15 



FIG. 8: Replotting experimental results of Ref. [TI| vs e. (a) 
Staying probability P. (b) A cff defined by Eq. ©. Fitting 
with Eq. © at e < 1 yields I ~ 11. 



going transitions at different values of the sweep field W 
become effectively decoupled. In all particular cases that 
we have considered, Iq can be estimated as the number 
of neighbors that are strongly interacting with a given 



TLS, \Jij\ > max(A, 



Vi - V, 



The sign of Iq depends 



on the sample shape in the case of the DDI. 

Our results differ from those of the mean-field ap- 
proximation. The MFA becomes valid for the descrip- 
tion of the LZ effect only in the case of weak interac- 
tions \Jij \ <C A/y^or in the case of superstrong disorder 
cry ^> max(| Jij |) if the interaction is strong. Our rigor- 
ous expression for Iq also describes models in which MFA 
predicts irrelevance of the interaction because of cancel- 
lation of molecular fields. Present results can be used as a 
"whetstone" for checking the quality of different approx- 
imations that can be suggested in the future. 

For long-range interactions Jij that exceed A over large 
distances, huge values of Iq are generated in the absence 
of the static disorder. In this case the region of fast 
sweeps, where a simple non-interacting LZ effect with 
P = 1 — e can be observed, becomes very narrow. The 
biggest values of Jo emerge in samples of ideal ellipsoidal 
shape, where the dipolar held is homogeneous. 

For samples of general shape gradients of the dipolar 



Acknowledgment 

The authors thank Wolfgang Wernsdorfer for supply- 
ing detailed information on the Fes crystal studied in Ref. 
[ill and E. M. Chudnovsky for useful discussions. 



APPENDIX: MEAN-FIELD APPROXIMATION 

Let us now consider the mean-field approximation for 
the many-body LZ effect at fast sweep and compare its 
results with our rigorous results obtained above. Here 
one has to solve the equations for single spins at different 
lattice sites in different effective fields: 



iMii = E tl a. a - —a,ii 

ih~a A = BifOif - y a a> (A- 1 ) 

where we used the definition of the wave function in the 
form of Eq. J3J at each lattice site i. Here the energies 
are 



1 

E% i — — 
2 



H z (t)-Vi + 2j2Jij(o-j 

j 

H z (t)-V l + 2j2-h (2M 2 - 1 



A.2) 



13 



and — —E^. Using the transformation similar to Eq. 
(|18|l one can rewrite these equations as 



da, 



dt 2h 



dt 



with 



where 



^ e -»i«i, 



*, (0) M 



vt 1 



(A.3) 



(A.4) 



(A.5) 



Vi is defined by Eq. (JTJJ, and 



(A.6) 



Now we can solve these equations for fast sweep e <C 1 
iteratively in A, similarly to Eq. I|20|l writing 



aa = 1 



2h °*i 



iA„(i) / zA\ _( 3) 
flit = «,t 



(A.7) 



Further one has to expand e l **W in Eq. I|A.3(I in |aj|(t)| 2 
that is small as s for fast sweep: 

e ««W = e ^ 0) (t) [i + + ...]. (A.8) 

Using Eq. (|31|) one can show that S&i contains the prod- 
uct £j9«. That is, making the expansion in powers of e, 
Eq. 10 within the MFA results in the expressions that are 
also automatically expanded in (3^. In contrast, within 
the rigorous formalism parameters e and /3y are split 
from each other and the terms of the e expansion con- 
tain general functions of (3^. Returning to the MFA, one 
obtains the expansion 



ag } (f) = 






J — OO 


- 




/ dt'e 


./ — OO 


a|?(*) = 


./ — OO 



a«(t") 



Now the staying probability at the site i 
P l (t) = l-|a iT (t)| 2 



(A.9) 



(A.10) 



can be expanded as 



2h 

A 

2h 



2 Re 



(a<?*(f)o<ffl),(A.ii) 



cf. Eq. I|22[l . In the limit t — > oo, the terms of this formula 
that do not contain Jy just reproduce the known terms 
of the expansion of the LZ probability e~ e in powers of 
e. One can rewrite it in the form 



Pi = l-e+ ^+e 2 E 7 ' 



(A.12) 



where 



J « = (^) 2Re 



and we have defined 



(a^M&gjM) (A.i3) 



6$(oo) = ^J« / dte-»5 0, W / dt' dW(t') 



'«TV^~ n^J^T Jo * T 

(A.14) 

To calculate Jy, it is convenient to introduce the di- 
mensionless parameters cti, fl^, and oy- defined by Eq. 
(I31|l and the dimensionless time variable 



2irh 



Then Eq. (|A. 5|> transforms to 

*?\t)=<p i (z)=n(z 3 -2a i z), 
and from Eq. 1JA.13J1 one obtains 

poo 

Iij = 47r/3 w Re /i (1) *(oo) / dze-^ z ) 



(A.15) 



(A.16) 



x i I dz' 
J o 



/fV) 



and 



jf } (z)= / dz'e^- (z,) . 

i/ CO 

One can integrate by parts in Eq. (|A.17|) : 

Iij = -4 7 r/3 y -Rc[z/ l (1) *(oo) 

x / dzf^(z) ff \z) 



(A.17) 



(A.18) 



(A.19) 



Here ft (z) can be expressed in terms of Fresnel integrals 
S(x) and C{x): 



dz'e-^*') = / dz'e-™( z ' 2 - 2S " z ') 
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X 

Z — Oii 



< ' I dz e 

-oo 



V2 



c 



iS 



V2(z - on) 

}■ 



With u = a/2 ( z — ay) one obtains, finally 



du 



(A.20) 



(A.21) 



C*(u - \/25ij) - S{u ~ y/Ydij) 



\ + C{u)\ 2 +(\ + S{u)~ 



For 8ij = this result simplifies to 



(A.22) 



*0i 
V2 J- 



du [C 2 (u) - S 2 (it)] = . ( A.23) 



Moreover, it can be shown that the symmetric part of 
Eq. i|A.22l) defined by Eq. I|25|l coincides with the small- 
Pij expansion of the rigorous quantum mechanical result, 
Eqs. © and (HO). 
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